Pinhole calculations of the Josephson effect in 3 He-B 
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We study theoretically the dc Josephson effect between two volumes of superfluid 3 He-B. We first 
discuss how the calculation of the current-phase relationships is divided into a mesoscopic and a 
macroscopic problem. We then analyze mass and spin currents and the symmetry of weak links. In 
quantitative calculations the weak link is assumed to be a pinhole, whose size is small in comparison 
to the coherence length. We derive a quasiclassical expression for the coupling energy of a pinhole, 
allowing also for scattering in the hole. Using a selfconsistent order parameter near a wall, we 
calculate the current-phase relationships in several cases. In the isotextural case, the current-phase 
relations are plotted assuming a constant spin-orbit texture. In the opposite anisotextural case the 
texture changes as a function of the phase difference. For that we have to consider the stiffness of 
the macroscopic texture, and we also calculate some surface interaction parameters. We analyze the 
experiments by Marchenkov et al. We find that the observed n states and bistability hardly can be 
explained with the isotextural pinhole model, but a good quantitative agreement is achieved with 
the anisotextural model. 



PACS: 67.57.De, 67.57.Fg, 67.57.Np 



I. INTRODUCTION 

The Josephson effects in superconducting weak links 
have been actively studied and applied since 1960's. By 
analogy, similar effects also exist between two volumes 
of superfluid connected by a weak link. There has been 
recent progress in observing the Josephson effect in su- 
perfluid 4 He |jj. In this paper we study superfluid 3 He, 
where the Josephson effect was experimentally confirmed 
over ten years ago by Avenel and Varoquaux [g,|| . How- 
ever, more recent experiments in Berkeley have raised 
new questions j|,[| which we are going to address here. 

Avenel and Varoquaux |2||3|] used a single narrow slit 
as the weak link in 3 He. They found current-phase rela- 
tions J{4>) that are very similar to those seen for tunnel- 
ing junctions or microbridges in s-wave superconductors. 
In these systems the relations are generally close to a sine 
function, J((f>) — J c sin (f>, or slightly tilted from this form. 
They are characterized by a single maximum supercur- 
rent J c and a negative derivative at <f> = 7r, J'(ir) < 0. 
The experiments at Berkeley used a 65 x 65 array of 
small apertures in 3 He-B ^j\^. At high temperatures 
J((j)) was found to be sinusoidal. At lower temperatures 
a "7T state" developed, where the derivative is positive at 
(j) ss 7r: J'(n) > 0. In addition, the weak link could be 
found in two distinct states with different current-phase 
relations. One of the "bistable" states had consistently 
higher critical current (H state) than the other (L state). 
Preliminary results of it states and multistability in a 
single narrow slit have also been reported || . 

Several theories have been put forward to explain these 
findings |7|-[ll||. It has been suggested that the reduc- 
tion of the Josephson coupling due to finite number of 
particles can lead to a it state in trapped atomic gases 



3 He were first calculated by Monien and 
. Their calculation used a very simplified 



p"2[ . This suggestion has been extended to 3 He [0,§. 
In the present paper we do calculations with the quasi- 
classical theory, which is an exact expansion in T c /Tp 
(the superfluid transition temperature over the Fermi 
temperature), and find no sign of such a mechanism in 
the leading order. Therefore we consider it unlikely that 
this mechanism could quantitatively explain the 7r state 
observed in 3 He. What look more promising are the- 
ories based on the 3x3 matrix structure of the order 
parameter in 3 Hc |iq , |ll^ , ^3|Jp| l . Unusual current-phase 
relations in 
Tewordt (l| 

one-dimensional Ginzburg-Landau model, and the physi- 
cal relevance of their intermediate branches around </> sa -k 
remains controversial. The first unambiguous evidence of 
a new branch in J(4>) came from 2-dimensional Ginzburg- 
Landau calculations Mm . Besides the usual case of par- 
allel n vectors on the two sides of the junction, this cal- 
culation considered also antiparallel n vectors, and the 
unusual J{<j>) was found only in the latter case. The new 
branch in J((f>) did not yet qualify as a i state, how- 
ever, because J'(tt) was found negative at the parameter 
values studied in Ref. [14|. More extensive Ginzburg- 
Landau studies in Ref. 1 1 ] found that a proper 7r state 
(J'(7r) > 0) occurs in the case of parallel n's through 
spontaneous symmetry breaking in a sufficiently large 
aperture. 

The Ginzburg-Landau calculations are applicable to 
relatively large apertures. A tractable opposite limit is a 
very small aperture, a pinhole. The pinhole model was 
first studied by Kulik and OmePyanchuk for an s-wave 
superconductor |15[ . At low temperatures J(</>) has con- 
siderable deviation from the sine form, but there is no 
7T state. Kurkijarvi considered the same problem in 3 He 



Jl6[ . In 3 He-B the order parameter is always modified 
near surfaces. Neglecting this complication, Kurkijarvi 
found that J(4>) for parallel fi's is exactly the same as 
for s-wave superconductors. Yip generalized Kurkijarvi's 
calculation to other orientations of the n vectors 10 1. He 



found a tt state for antiparallel fi's, as well as for some 
more complicated configurations, which can occur in a 
magnetic field H > 1 mT. We call this mechanism of the 
tt state isotextural because the texture [the field n(r)] is 
kept constant while calculating J {(f). 

The discussion above concerned a single aperture. 
There exist three different suggestions as to how a tt 
state can appear in an array of apertures. Avenel et 
al. assumed that if the individual apertures have a hys- 
teretic J(<j>), approximately half of the apertures could 
be on a different branch than the others. The net effect 
would be the formation of a tt state ||. We consider 
this explanation unlikely because apparently the aper- 
tures in Ref. M are not hysteretic, and also because it 
is difficult to understand why exactly half of the aper- 
tures could behave differently. The second alternative is 
that the tt state appears trivially in a coherent array if 
an isotextural tt state appears in each of the apertures 
independently. The third alternative is an anisotextural 
tt state, where the n texture changes as a function of <j) 
p| . This mechanism can lead to a tt state even in the 
case when the isotextural J(<f>) is sinusoidal. 

The purpose of this paper is to study J{<j>) in 3 Hc-B 
as completely as possible using the pinhole model. Sec- 
tion pi starts with a division of the problem into meso- 
scopic and macroscopic, and using symmetry arguments 
derives general expressions for the Josephson coupling. 
The mesoscopic problem is discussed in the following four 



sections. Sec. [II introduces the quasiclassical theory and 



the assumptions relevant for our calculation. The pinhole 
model is defined in Sec. fV\ a nd a general pinhole energy 
functional is derived Sec. |V]. The propagators are cal- 
culated in Sec. 



VI 



and the Josephson energy and the 
currents are evaluated. This corrects the calculations by 
Kurkijarvi |L6| and by Yip [[10| by using a selfconsistently 
calculated order parameter. We consider both diffuse and 
specular reflection of quasiparticles at the wall. We also 
discuss the case where scattering is present within the 
pinhole. 

The rest of th e pa per is devoted to the macroscopic 
problem. In Sec. VII we discuss the interactions that are 
important on the macroscopic scale. We estimate the 
length scales and evaluate some surface-interaction pa- 
rameters that have not been calculated before. In Sec:. 



VIII we plot isotextural current-phase relations in dif- 
ferent situations. The anisotextural Josephson effect is 
discussed in Sec. [X. The effect is first demonstrated with 



a simple model. Then we estimate the textural rigidity, 
and calculate current-phase relationships for an array of 
pinholes. Section DO is devoted to the analysis of different 
tt state models in the Berkeley experiment H. Section 




FIG. 1. A weak link between two bulk superfluids, L and 
R. The arrows denote the n vector of 3 He-B. The dashed 
line divides the liquid into a mesoscopic region (inside) and a 
macroscopic region (outside). 



XI finishes with some conclusions. 



II. SYMMETRY CONSIDERATIONS 

The Cooper pairs in superfluid 3 He are in a relative p 
wave state and form a spin triplet. This state of affairs is 
reflected by the 3x3 tensor character of the order param- 
eter A^. The first index in A^i (greek letter) refers to 
the spin states and the second index (latin letter) to the 
orbital states. In the bulk of 3 He-B the order parameter 
has the form fl7|] 



A^i = i? M iAexp(k/ 



(1) 



Here A is a real constant and R^i is a rotation matrix 
satisfying R^R^j = Sij. (We shall consistently assume a 
summation over x, y, and z for repeated index variables 
/i, v, i, j, a, /3, 7, etc.) The rotation-matrix structure 
and A are fixed on the scale of the superfluid condensa- 
tion energy. On this energy scale, the state (m remains 
degenerate with respect to the phase cf) and different ro- 
tations R^i, which are parametrized with an axis n and 
an angle 8. The degeneracy with respect to <f> and R^i 
(or equivalently <p, n, and 6) is partly lifted by weaker 



interactions, which are discussed in detail in Sec. VII. 

Imagine now two bulk volumes of 3 He-B connected by 
a weak link, see Fig. [y. We consider the system to consist 
of a mesoscopic region at the junction and a macroscopic 
region outside. There can be considerable variation of the 
soft degrees of freedom in the macroscopic region, as illus- 
trated for the n vector in Fig. [|. The mesoscopic region 
is chosen sufficiently small so that (i) both <p and R^ are 
effectively constants at the macroscopic-mesoscopic bor- 
der, and (ii) the weak interactions affecting <p and R^i 
can be neglected in the mesoscopic region. The size of 
the mesoscopic region is limited from below by the con- 
dition that the bulk form of the order parameter (||) has 
to be valid in the macroscopic region. Thus the size of 
the mesoscopic region has to be large in comparison to 



both the superfluid coherence length £o ~ 10 nm and the 
size of the aperture. The mesoscopic region can also be 
chosen to cover several apertures. 

The rotation matrices and the phases on the left (L) 
and right (R) sides are generally different. Their values at 
the macroscopic-mesoscopic border are denoted by R \ 
and <f> L ' . The most general form for the Josephson cou- 
pling energy associated with the weak link is given by 



F, = F, 



b R ,R L 



R 

fit ■> -^uj 



(2) 



However, due to global phase invariance there can really 
be only a dependence on the phase difference cp = (p R — 
cf> L . In addition, if we assume the intervening wall to 
be magnetically inactive, then we should have invariance 
with respect to global spin rotations as well. This means 
that the energy can only depend on the product of the 
rotation matrices, or the quantities ipij — R^Rfy, so 
that 

Fj = F 3 (cfr,R%RX). (3) 

Using the functional Fj (0) one can now calculate two 
conserved currents. Firstly, the mass current through 
the aperture is given by 



J 



2m 3 dFj 



(4) 



we fix the z coordinate to be along the axis of the weak 
link. Now, if the twofold rotation symmetry around z is 
replaced by a fourfold symmetry ( — — — ), then an ex- 
change of ip xx and ipyy must not affect Fj . Finally, if the 
rotation symmetry around z is continuous (^^), the 
dependence can only be through ip zz and the invariant 



combination ip x 
Fj = Fj( 



+ ipyy, that is 



'mi 

r>L r>B 
• -^ fix -^ fix 



^liy fJ,y ' n 'fj,z Ix fiz ) • 



(6) 



Close to T c the amplitude A of the order parameter 
(|l|) approaches zero, so that we can expand Fj(A^ i , A^j) 
in powers of A. In order to be consistent with Eq. ra), 
the leading order term in the expansion must be 

Fj = -[aRfutfL + 0{R^ x R% + tfvtfv)] cos <P, (7) 

where a and j3 are some real valued phenomenological 
constants. In Ref. fiTj this was introduced as the Joseph- 
son energy of the tunneling model, but as the derivation 
above shows, it is more general. The tunneling barrier 
for 3 He was first considered in Ref. [f20| . 

In order to further determine the functional Fj (||,f||3) 
it is necessary to do a calculation in the mesoscopic re- 
gion. This is discussed in the following sections |t|JVI|. 



III. QUASICLASSICAL THEORY 



See, for example, Ref. |fj| . Secondly, due to the different 
spin-orbit rotation matrices on the two sides, there is also 
a spin current flowing through the aperture. Likewise, 
this is obtained from Fj by differentiation: 



rspin _ rth jdR 



OF, 



fii WXY 



(5) 



This can be seen by replacing ipij in (||) by R^R^Rap 



with the relative rotation R a n = 5 a ff+e 



e a/ 3 7 0f 7 , 



and iden- 



tifying SFj = J*P in S6 y + 0(56^) |9[. 

The energy and the currents also satisfy some other 
symmetry relations. For example, assuming invariance 
with respect to time reversal, we have i ?1 j(— <fi, ipij) = 
Fj(<p,ipij). As a consequence J{—<p,ipij) = —J(<p,ipij) 
and J^ in (-0,7/.y) = J*v ln {<p,ipij). Also, since the 
phase factor exp(i^) is defined only modulo 2-7T, we have 
Fj(4> + 2ir, ipij) = Fj((p, ipij). Here we must keep in mind 
that in a long junction Fj is not a single- valued function. 
Analogously, there is periodicity with respect to the rota- 
tion angle, which is automatically contained in the form 
of Eq. (g), since i?^i(n, 6 + 2tt) = i? M i(n, 6). 

The functional form of Fj (g) can be restricted further 
if we assume some additional symmetries. For example, 
if the aperture is symmetric under a parity operation, 
we have Fj(<p, ipji) = Fj(<p,ipij) and thus J* pin (<p,ipji) = 



. ppi: 



, ipij ) . If the aperture has full "orthorhombic" 
symmetry — — — , then Fj can only depend on the rota- 
tion matrices through ip xx , ip yy , and ip zz . From here on 



We use the quasiclassical theory of Fermi liquids |21| to 
calculate Fj (||) for a pinhole aperture. Here we present 
the theory only in the depth needed for the following. 

The central quantity is the quasiclassical propagator 
g. In the stationary case which we are considering this 
can be written as i?(k, r, e m ), where r denotes spatial po- 
sition, k parametrizes a position on the Fermi surface, 
and e m — ■Kk^T(2m + I) are the Matsubara energies. 
The propagator is determined by the Eilenberger equa- 
tion and the normalization condition 



[ie m T 3 - (7, g] + ihv F k ■ V r g = 

99 = -1- 



(8) 
(9) 



Equation (||) can be interpreted as describing transport of 
quasiparticle wave packets which travel on classical tra- 
jectories with the Fermi velocity vp = i>Fk. The propa- 
gator g as well as the self-energy <y(k, r) are 4x4 matrices, 
reflecting the spin and particle-hole degrees of freedom of 
a quasiparticle; the matrices r% (i — 1, 2, 3) are the Pauli 
matrices in the particle-hole space. 

The 2x2 spin blocks of g can be decomposed into 
scalar and vector components as 



9 



5 + g;5i (/ + f • 0L)i£ 2 
i£ 2 (/ + f - £) 9~ 3L2S ■ °L2.2 



(10) 



where g_ = xa^ + ya 2 + zer 3 , and a_ t (i = 1, 2, 3) are the 
spin-space Pauli matrices. The self-energy a is written 
similarly 



v + v ■ a_ A • qyo_2 
ig_ 2 A* • er v — a^u ■ a c, 



■ u_ v — v_2 ° 



(11) 



Here the off-diagonal terms contain the p-wave pairing 
interaction in the form of the gap vector A M (k, r) = 
A i 
equation 



nt( r )k%- This is determined by the self-consistency 



rk B Tj2 



* 3 /*f(kV, em )(k'.k) 



T 
-A In — = 0. 



(12) 



where T c is the superfiuid transition temperature. This 
form is valid in the weak coupling approximation, where 
the quasiparticle-quasiparticle scattering is neglected. 
The diagonal components v, v and f, i> are the "molec- 
ular field" self-energies arising from redistributions of 
quasiparticles. The scalar parameters v and v arise in re- 
sponse to mass currents, and they turn out to be negligi- 
ble as will be argued in Sec. |v|. The (real valued) vector 
parameters v and u describe the response to a magnetic 
field or spin currents. As discussed above, the magnetic 
field can be neglected in the mesoscopic region. In con- 
trast, there are always rather strong spin currents flowing 
along surfaces in 3 He-B |22J. These have to be taken into 
account with the self-consistency relation (£> = — v) 



u(k,r) = nk B T 



Z^ 



dfifc, 

47T 



A Q (k-k')g(k',r,e m ). (13) 



Here A a (x) = E^o^ t 1 + F l '/( 2l + l)] _1 fl(aO> Pi are 
the Legendre polynomials, and all terms with even I drop 
out due to symmetries. 

For our purposes the most important physical quanti- 
ties to be evaluated from g are the mass supercurrent 

j(r) = 2m 3 v F N(0)7rk B Tj2j ~jf k 9( k > r, e m ) (14) 

m 

and the spin supercurrent 
j^ pm (r) = Hv F N(0)irk B Tj2 / -^k ff7 (k,r,e m ). (15) 

m ■* 

Here N(Q) = m* 2 vp /(2ir 2 h 3 ) is the one-spin normal- 
state density of states at the Fermi surface, where m* — 
1713(1+ F*/ 3) is the effective quasiparticle mass, 1713 being 
the mass of a bare 3 He atom. The superfiuid coherence 
length is defined by £0 = fiv-p / (2irk B T c ) . For Ff, vp, T c 
and other pressure dependent quantities we use the vapor 
pressure values whenever needed. In Eq. ( |13[ ) we assume 
Ff 1 = for all odd I > 3. Since the parameter Ff is not 
well known, we usually set it to zero also, but values in 
the range — 1...0 have been tested. 




AR? -J* 



FIG. 2. The mesoscopic region (Fig. hi) for a pinhole. Two 
quasiparticle trajectories are shown. The coordinate u is plot- 
ted along the straight transmitting trajectory. The arch on 
the left hand side of the pinhole denotes an imaginary surface 
that is used to close the pinhole (Sec. M). The diameter D 
of the pinhole and the thickness of the wall W are assumed 
small in comparison to £o- 



IV. THE PINHOLE PROBLEM 

Consider the case of a single pinhole in a thin wall 
separating two volumes of 3 He-B - - the situation of 
Fig. |[ The hole can be thought of as a pinhole (i.e. 
"very small") and still be treated quasiclassically, if its 
dimensions (diameter D and wall thickness W) satisfy 
£0 ^> D, W 3> Af, where Af is the Fermi wavelength. In 
addition 3 He is in the pure limit, where the mean free 
path of quasiparticles I ^> £o- Usually one further as- 
sumes that W -C D so that scattering in the aperture it- 
self can be neglected. We allow for a finite W/D and thus 
consider also deflected trajectories of the form shown in 
Fig. @. The pinhole limit was first considered within the 
quasiclassical theory by Kulik and OmePyanchuk in the 
case of s-wave superconducting microbridges Eq ]. Sev- 
eral previous calculations for the spin-triplet case of 3 He 
also exist |bj,[ll],[16 23 26| . For another type of quasiclas- 



sical Josephson model, see Ref. [|27j . 

What makes the pinhole case attractive is that no self- 
consistent calculation of the order parameter in the aper- 
ture is needed. More precisely, the leading term in the 
coupling energy Fj is on the order of the superfiuid con- 
densation energy in the volume D 2 £q. This effective vol- 
ume is large compared to the volume (~ D 3 or WD 2 ) of 
the pinhole. The leading term in Fj (and thus also the 
corresponding terms in the currents) can be computed 
using a (O) that is calculated for a planar wall without 
the pinhole. There is a correction ex D 2 to a, but because 
of the stationarity of the energy functional with respect 
to a, this affects Fj only in the order oc D 4 /£o, which is 
negligible for a pinhole. Therefore we can use the order 
parameter profiles calculated for a planar wall. 

Determining the suppression of the order parameter A 
at the wall is thus the first step needed for our calculation. 



In the absence of mass currents and magnetic scattering 
it is sufficient to consider the parametrization 

A(°)(k,z) = [A ± (z)zz + A||(z)(xx + yy)].k, (16) 

where z is perpendicular to the wall. The gap functions 
A±(z) and An (z), which are real valued, are calculated 
self-consistently as explained in Ref. |E3|. For g we use 
the "randomly oriented mirror" (ROM) boundary condi- 
tion at a specular or a diffusive surface |2q]. The numer- 
ical calculation of g is described in Sect. [IV|. Our results 
for Aj_(z) and A|| (z) are similar to those found previously 
with ROM and other surface models [g2|,£5],£6 28 3_jj . In 
order to incorporate different phases and spin-orbit rota- 
tions on the two sides we write 



A(k,z) 



exp(i<j> L )R L ■ A(°)(k, z) for z < 
exp(i^ fl )^- A(°)(k, z) for z > 



(17) 



The thin wall is located at z — and A' ) §w£) is assumed 
to be symmetric: A±(—z) = A±(z) and Aii(— z) = 

A|,(z). 



V. ENERGY FUNCTIONAL 

Here we derive a general quasiclassical expression for 
the Josephson coupling energy in a pinhole. The deriva- 
tion follows closely the lines of a quasiclassical treat- 
ment of impurities in 3 He or superconductors p4p3[ . We 
start from an expression for the energy difference between 
states with one impurity and no impurity, V being the 
impurity potential |2l| , |34[| . For a small spatial range of 
V, the self-energy E can be assumed to be unchanged by 
it, and we get |32|] 

Sfl tot = -~ Tr[ln(-(V + S + V) - \n{-G^ 1 + £)]. 

(18) 
The trace operation Tr is defined as 

x-^ f d 3 k 

TrF = fc B T^ / —-Tr 4 ^(k,k,e m ), (19) 

where Tr 4 denotes the trace of the 4x4 Nambu matrix 
F . To eliminate the logarithm, we may apply some form 
of the "A-trick" |3q] . We choose to integrate over the 
strength of V by making the substitution V —f XV and 
writing 



Stf 



l Tl J ^(G^-t-XV)-'XV 



^Tr r^f, 



(20) 



The latter equality follows from a formal application of 
the t- matrix equation T\ — XV + T\G\XV and the rela- 



tion G = Gi + GiT x Gi. Here G = (G|7 



XV)- 1 



gives the full propagator in the presence of an impurity 
scattering potential. The intermediate Green's function 
G\ corresponds to a quasiclassical g± which has no dis- 
continuity at the impurity. Equation (EOT) is now in a 
form where the propagator can be |k|-integrated directly. 
However, to avoid a divergence in the Matsubara sum- 
mation, we have to subtract from (Bfl) the normal-state 
contribution i5fi N , obtained by setting E = E N in (p"8]). 
We define <5f2 = <5f2 tot — 50 N and transform this to the 
quasiclassical form 



SSI = ^hv F N(0)nk B Tj2 / -47 / 



x dA 

A 



xTr 4 [gi(k,ri m p,e m )i A (k,k, e m ) 

~9l (^1 r imp, £m)t\ (k, k, e m )J, 



(21) 



where ri mp is the location of the impurity and £.\(k, k, e m ) 
is the forward-scattering part of the t matrix, which is 
obtained by solving 



t A (k, k', e m ) = Aw(k, k') + Att^(O) 



*£&* 



x<?i(k", r imp , e m )t A (k", k', e ri 



(22) 



The energy formula ( pl| ) is simpler to use than those in 
Refs. |32||33|], since it does not involve an integration over 
r. More importantly, g\ is constant in the A integration. 

Now we specialize the above approach to the pinhole 
problem. The coupling energy Fj we wish to know is, 
by definition, the difference in energies between an open 
pinhole and a blocked pinhole. It should be irrelevant 
how the hole is blocked as long as the transmission of 
quasiparticles is prevented: changing the type of blockage 
should only change some constant terms in the energy, 
which do not depend on the soft degrees of freedom in (||). 
We choose to block the pinhole by a surface just on the 
left hand side of it (Fig. 0) . This surface is now considered 
as the "impurity" in Eqs. (|l8||2^). The coupling energy 
is then equal to Sfl ( |2l| ) except for a minus sign: _Fj = 
— dfl. The intermediate g\ is the exact propagator that is 
calculated for an open pinhole, and we drop the subindex 
1 from here on. 

The simplest choice for the blocking wall is a specularly 
scattering surface. It corresponds to a delta-function 
scattering potential V<5(1 • r) for the (infinitesimal piece 
of) flat surface 1 • r = in the limit V — ► 00. The t ma- 
trix of this type of impurity is of particularly simple form 

MM 



i A (k,k',e rn ) = 



2v F |k-l|AVdA5. 



(2) 



2w F |k-l| - AV[g(k, r imp ,e m ) + <?(k,r im p, e m )] ' 



(23) 



where dA is an area element (A|* <§C dA <C £q) of the 
blocking piece of wall with normal 1. The component of 
k parallel to this wall is denoted by k|| = k — (k • 1)1, 
and k = k — 2(k • 1)1 is the mirror-reflected direction. Wc 
insert i\ (|23| ) into Sfl (|2l| ) and integrate over the blocking 
surface. Performing also the A-integration, we find 



*j(V) 



x In 



1 



2 

2v F |k-l 



d^E 



^ F -/V(0)7rfe B T 

V[g(k, r imp ,e m 



dfi£ 

47T 



Ik • l|Tr 4 



.9\— 7 I"i m p, € m JJ 



2 UF |k-l| + 2iVr 3 Sgn(e m ) 



(24) 



where the normal-state propagator g N — 1T3 Sgn(e m ) was 
used. Next we take the limit V — > 00 and use the 
general properties Ti^lnA = lnDet 4 A and Det 4 A.B — 
Det4^4Det4-B, noting further that Det 4 [iT3 Sgn(e m )] = 1. 
All trajectories that do not pass through the hole in the 
absence of the blockade can be neglected, and thus the 
integration over the surface A can be transformed to one 
over the cross section (of area A ) of the hole. This leads 
to the general pinhole energy functional 
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where k z = k • z, the brackets (...) denote average over 
the trajectories (at fixed k) that hit the area A of the 
pinhole (A — irD 2 /4 for a circular hole), and L denotes 
the location immediately on the left hand side of the hole. 

The coupling energy (^5|) depends on the shape of the 
blocking wall only through the directions of the reflected 
momenta k = k — 2(k • 1)1. Because the surface can 
be chosen in different ways, there is a lot of freedom in 
choosing the reflected momenta. A particularly simple 
choice is retroreflection: k = — k. This can be achieved, 
for example, by a scmispherical blocking surface of ra- 
dius R satisfying £ ^> R ^> D, centered at the pinhole. 
This choice can simplify practical calculations consider- 
ably, since one can use symmetries effectively. After this 
one may, for example, calculate the determinant by using 
Det 4 (gi + .g 2 ) = [Det 4 (-21 + {gi,h})} 1/2 - 

The choice to block the surface on the left hand side 
was arbitrary and the blocking could equally well be done 
on the right hand side. The reason for not blocking in the 
middle of the hole is that there may be scattering taking 
place there. For example, a quasiparticle may hit the wall 
inside of the hole and be deflected (Fig. 0). In general 
this causes g to be discontinuous at the hole, and in order 
to give an unambiguous prescription for Fj (J25|), one has 
to specify one of the sides. If there is no scattering in the 
pinhole (W/D = 0), the propagator is continuous and 
the energy functional can be evaluated in the middle of 
the hole. In the absence of scattering also the average 
(. . .) in Fj (25) is trivial and can be dropped. 



Apart from the assumption of a pinhole aperture, the 
energy functional (Ea) is very general. Below we shall 
apply it for the special case of 3 He-B. 



VI. PINHOLE CALCULATION 

A. Trajectories 

For the case W/D = and no scatterers localized 
within the pinhole, all trajectories hitting the orifice are 
directly transmitted. For the case of finite W/D we con- 
sider a model which is based on the ROM boundary con- 
dition n). 

The ROM model assumes that the surface consists of 
microscopic pieces of randomly oriented mirrors. There- 
fore, any trajectory hitting the surface is simply deflected 
into another direction and the physical propagator is con- 
tinuous along it. When this process is averaged over a 
length scale which is large compared to the size of the 
mirrors (Af <C mirror size <C W, D), one obtains a prob- 
ability distribution for the scattering from one direction 
to another. Consider a quasiparticle coming out of the 
pinhole in direction k. The probability density that it 
entered the hole from direction k' can be written as 



Wtfr =p(k)fc t , + (l-i»(k)W t ,. 



(26) 



Here p(k) is the probability for direct transmission, and 
w k i.' ^ s * ne scattering distribution obtained by averaging 
over the surface of the pinhole which is visible from direc- 
tion k. For given k the distributions are normalized ac- 
cording to /(dfifc,/47r)Wk£, = 1 and f(dQk//47r)«?kk' = 
1, where the integrals are over all directions (backscatter- 
ing is also possible) . For a circularly cylindrical aperture 
of diameter D in a wall of thickness W one finds a simple 
form for the transmission function p(k) : 



pW 



^(7 — cos 7 sin 7) for ■& < arctan(Z?/IF) 
for ■& > a,icta,n(D/W), 



(27) 



where 7 = arccos(VT/D) tan $ and $ is the polar angle of 
the trajectory (k • z = cos$). 

Calculating the w^ ?. for a general type of ROM 
surface and a finite ratio W/D is difficult, since one 
has to consider multiple scattering. We shall restrict 
to a cylindrical aperture and start with fully diffu- 
sive, rough walls. In this case one can in principle 
expand w^, = P(k|k') = / d 2 riP ou t(k|ri)P in (ri|k') 

+ / d 2 n J diriPoutik^Pirilr^Pininlk') + • ■ •, where 
the Ti parametrize positions on the surface of the cylin- 
der. For given outcoming direction k the function 
Pout(k|r) gives the distribution in r for its origin, whereas 
the function Pj n (r|k') gives the distribution in incoming 



k', given a position of impact r. The diffuse interme- 
diate scatterings follow a "Markovian" process, so that 
the distribution P(r2\r\) is independent of the incoming 
direction. The formal expansion parameter here is W/D 
and the first term in w k k , is of zeroth order in it, corre- g 
sponding to a single scattering event. The second term 
has one intermediate scattering and is thus of first order, 
and so on. In the limit W/D — > only the zeroth order 
term remains, and the functions P[ n and Pout are assumed 
to approach simple "cosine laws": P m (r\k') oc |s(r) • k | 
for s(r) • k' < and P out (k|r) oc |k • s(r)| for k ■ s(r) > 0, 
where s(r) is the surface normal at position r. Upon nor- 
malization and insertion into the expansion for w^ k , one 
finds 



w k,k' = sin,9 k'( sin Xk,k' 



x k ,k' C0S Xk,k'): 



(28) 



where x k k' ^ s ^ ne difference in incoming and outgoing 
azimuthal angles and $ k , is the incoming polar angle. 
This distribution is largest for angles near Xfc k' = n an< ^ 
i9 k , = 7r/2, i.e., for scattering into directions close to the 
plane of the surface. 

As a second case, consider the possibility of specular 
scattering in the pinhole. Then a k' directional quasipar- 
ticle hitting the surface at position r will reflect into the 
direction k = k — 2(k' • s(r))s(r). In this case the pre- 
vious distribution functions P out , Phi and P have to be 
generalized a bit to take into account the non-Markovian 
character of the scattering, but in the limit W/D — > 
no problems will arise. Since s is in the xy plane, we 
necessarily have k • z = k' • z, and a similar calculation 
as for the diffusive case gives 

w kk , = TrCsin^)- 1 ^ - tf k ,) sin( Xkik ,/2). (29) 

More refined distributions could be obtained by taking 
into account higher-order terms in the expansion, but 
doing so analytically would be difficult. In Ref. |27J these 
were briefly discussed in the case of a long pore with 
specular walls. 



B. Propagator 

Here we describe briefly the method used to generate 
the propagators. The idea is to calculate the propagators 
numerically only for A^ ^ (16), and to obtain analyti- 
cally the dependence of the true propagator on the soft 
degrees of freedom in A (|7|). The matching of the left 
and right solutions at the pinhole is most conveniently 
done with the "multiplication trick" ]22| , p5| , p3[ | . There 
one first calculates two unphysical solutions g< and g> 
of the Eilenberger equation (ph and g is constructed using 



Here we denote by g < and g> the solutions decaying ex- 
ponentially towards left (z — — oo) and right (z = +oo), 
respectively, independently of the direction of k. 

We rewrite the propagator components as g = c + d, 
= c + d, f = a + b, f = a + b and g — c — d, g = c d, 
/ = a — b, f = a — b. In terms of these, the Eilenberger 
equation decouples conveniently into three independent 
blocks of linear, first-order differential equations which 
are numerically more convenient to handle ]22| . The first 
task is to find the unphysical solutions consisting of com- 
ponents a, 6, and c. For the real valued order parameter 
A' ', the unphysical propagator components can be cho- 
sen such that a and b are real and c is purely imaginary. 
The unphysical block of equations p2| then becomes 

£ m b + -hv F d u a = 

e m a + A (0) • Imc + -hv F d u b = 

v x Imc + A (0) 6+ -hv F d u Imc = 0. 



(31) 
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g(i.,r,e m )=iS^) ^ r ^^ r ^ . (30) 
{ff<(k,r,e m ),5>(k,r,e m )} 



Here u is the coordinate along an arbitrary trajectory 
r = ro + uk, and we fix u — at the wall (z — 0). In 
accordance with (|3(j), the exponential solutions of ( |3l| ) 
which go through the pinhole are denoted by <7<, and 
gy , respectively. These are the solutions that are natu- 
rally obtained by integrating from the bulk towards the 
wall on L and R sides. Because of symmetries we only 
need to calculate g K , and we do this using fourth-order 
Runge-Kutta method. We introduce a short-hand nota- 
tion for the numerically calculated quantities 



(32) 



These were evaluated for several directions of k, whose 
polar angles ■& were chosen so that angular integrations 
in ( p8[ ) and ([Io|) could be carried out using the Gaus- 
sian quadrature, usually with 32 points in the range 
cos?? = — 1...1. (Due to the symmetry of the inte- 
grands, only values for k z > actually need to be con- 
sidered.) The number of (positive) Matsubara energies 
e m = irk B T(2m + 1) was between 10 to around 100 de- 
pending on the temperature. Above we assumed A'°) 
and v to be already known. The method for their self- 
consistent calculation is the same as above, except that 
also the solutions diverging away from the wall have to be 
calculated. The initial conditions for these were obtained 
using the specular reflection or the ROM boundary con- 
dition. 

The functions <?> can be obtained by using the rela- 
tions 



aX,' (k, u, e m ) = +a ( °> (k, -u, e m ) 



bxP (k, u, e m ) = -b { °> (k, -it, e m ) 

c^ 0) (k, u, e m ) = +c^ 0) (k, -u, e m ), (33) 

which are based on the symmetry 

A (0) (k, M ) = A (0) (k,-u). (34) 

From the solutions g<- for A*- ' we get the solutions g < 
for the general A ( JL7| ) on the L side by forming the linear 
combinations 
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(35) 



The same equations hold on the R side when L is replaced 
by R and < by >. 

The physical propagator at the pinhole (r = ro = 0) 
can now be obtained using Eqs. (p0|), (|32|), (|33"|), and (35). 



For the case of deflected trajectories we have to specify 
separately the momentum k' on L side and k on R side. 
Only the transmitted trajectories (k z k' z > 0) need to be 
considered, and we get 

a k,k' ( £ ™) = i S k.L t C ' L ^ lA ain 2^^ B C0S 2 ^ 
+(L4' sin -<f) - S' cos -0)C fl ] 

b k,k'( e ") =is k,y c i (^cos-0-LBsin-0) 



-(A' cob -^-iB' sin -0)0*] 



dtt-(er 



i s. 1 . [i(AA' + BB 1 ) sin - (45' + A'B) cos < 



d k,k'( e ") =-*tLC£ x C R , 



(36) 



k,k' 

where Cl,r = R lr - C, and primes denote values corre- 
sponding to direction k'. The normalization constant is 
given by 

s k,k'( £ ™) = ~ s -k,-k'( e ™) = s & n C k z) (37) 

x [-{AA' + BB') cos <t> + \{A'B + AB') sin cj) + C' L - C R ] . 

For the case of direct transmission (k 7 = k) these expres- 
sions simplify considerably. 

C. Currents and coupling energy 

As an application of the above results, consider the 
Josephson current in the pinhole. Using the general sym- 
metries for propagators, the mass current density (Tbl) 
can be written in terms of Red alone. The total cur- 
rent is then given by Re d k k , (£56|) integrated over the 
distribution (J2q) of trajectories: 



J = A o 2m 3 v F N(0)nk B T 

'dfit, 



E/SM 



j„ , 4l -»t,k.R<" i k,k-(^). (38) 



In the case of direct transmission only (W/D — or 
W k k , = (5 k k , ) , one can apply trigonometric identities to 

put d(k, 0, e m ) = d k k (e m ) in the form 



d(k,0,e m ) = Sgn(A;,) 

1 ^ (S 2 -A 2 )sin(0 + aC) 



4E 



2UB 



4^ 1 ^ 2 sin 2 [i(0 + aC)]+i? 2 cos2[i(^ + aC^ 



(39) 



where C(k, e m ) is defined by Cl • Cr = C 2 cos£. The 
real part of ( |39| ) is now equivalent to equation (1) in Ref. 
Eof , but more general. The quantities 0±£(k, e m ) are the 
effective phase differences experienced by quasiparticlcs 
with different spin projections along the axis C^ x C R 
oc Aj,(k) x A^(k). In the special case that A is assumed 
constant, i.e., unsuppressed at the walls, the Matsubara 
summation can be done analytically and one obtains the 
same result for J as in Ref. fllCf ]. 

For the Josephson coupling energy (|25j) we find 

Fj = ^A c ho F N(0yxkBTj2 f~^\k\ (40) 

m 

^^k.k'l^l^k'^OP-lnWAA') 2 ]}. 

In the first term there appears the squared modulus of the 
trajectory-invariant normalization constant (pTJ). The 
second term has to be retained to have convergence in 
the Matsubara summation. 

Consider again the direct-transmission case W/D = 0. 
In the Ginzburg-Landau limit T — > T c we can verify the 
phenomenological form (0) and calculate the parameters 
a and (3. In this limit the amplitude of A is small 
and since |A| 2 ~ |C| 2 = \A 2 — B 2 \, we should have 
\A 2 - B 2 \ < A 2 + B 2 , A 2 » B 2 and |s kk | 2 « AA 4 . 
It follows that the logarithm in the first term of (|4(j) can 
be expanded to linear order to give (R), where 
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(41) 



Figures shows the temperature dependence of the tun- 
neling parameters for a diffusive wall with A chosen to 
match the total open area of a coherent array of holes 
with dimensions as in Ref. H . Also shown is the textural 
rigidity parameter 7 whose role is to be discussed below 
in Sect. IX. Close to T c the strength of the coupling, 
i.e., the parameters a and j3 go as a, f3 oc (1 — T/T c ) 2 , 
whereas the rigidity 7 ex (1 — T/T c ). 
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FIG. 3. Temperature dependence of the tunneling model 
parameters a, (3 (Ww, and 7 (B7I) as calculated for a diffusive 
wall and total open area A — Ak where, A — 3.8 • 10 -8 m 2 , 
k = 14.7 ■ 10~ . The pressure is bar. 
for i — 1, 3, etc. 



W/D = 0, and Ff = 



The Josephson current ( |3§| ) and energy (40) were ob- 
tained completely independently It is essential to check 
that they are consistent with each other. One can easily 
see that the component Re d^ ?, satisfies 

Red^,(e m ) = Isgn(fc z )^ln| Sfcic/ (e m )| 2 , (42) 
and hence the macroscopic current formula (0) is exactly 



satisfied. As a further check of the energy (40) we can see 
that also the spin current formula (ra) is satisfied. Using 
C^ • C fl = R^RKjC-Cj the energy @ is seen to be a 
function of the products R^R^ and calculating the spin 
current from (p|) is thus possible. Writing also [C L x 
C_r] 7 = e a p 1 R^ k R% l C' k Ci in the propagator component 
(36) and using the quasiclassical spin current expression 
(15) it can be checked that the result for Jf; pin agrees 
with the one obtained from (|5[) . 



A. Interactions and coupling constants 

The most important hydrostatic energy term arises 
from the dipole-dipole interaction between the nuclear 
magnetic moments, 



F D = 8 5D A ; 
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where 9 is the rotation angle of the spin-orbit rotation 
R^i(h, 8). The effect of Fe> is to fix 9 to the value #0 ~ 
104° in the bulk liquid. 

There is no conflict between Fd ( [43| ) and Fj ^) . Even 
if R I have their rotation angles fixed to 9q, Fj depends 
on their product ipij — also a rotation matrix — and it 
can attain all possible values if n L and n R are directed 
properly. Thus 9 is not changed from 9q by the Josephson 
coupling. The same applies to all surface energies below. 
We therefore assume Fr> to be in its minimum everywhere 
and study only the position-dependence of n(r), which is 
known as the texture. 

In the absence of a magnetic field the dominant inter- 
action determining the texture near a wall is the surface- 
dipole interaction 



F S b= /d 2 r[6 4 (s-n) 4 ~6 2 (s.n) 2 ]. 
Js 



(44) 



where 62 and 64 are positive coupling constants and s is 
the surface normal. There are usually many walls with 
different orientations present and therefore there is a con- 
flict between their orienting effects. This leads to gradi- 
ent (bending) energy 
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(45) 



which is related to spin currents in the bulk. The gradient 
energy also has a surface part 



■Fsg = A S g / d 2 rsjR aj — — 



(46) 



VII. TEXTURAL INTERACTIONS 

As discussed in Sect. O, the Josephson effect in 3 He 

— T R 

depends on the rotation matrices R ■ on the two sides 
of the weak link. These matrices are determined by the 
competition of a number of relatively weak bulk and sur- 
face interactions, which lift the degeneracy of the B phase 
order parameter (|lj). The equilibrium configuration is 
found by minimizing a hydrostatic energy functional. We 
shall present the hydrostatic theory to the extent needed 
here. For a recent review see Ref. |36|1. 



In the presence of a magnetic field another surface inter- 
action becomes important: 



(47) 



F S n = -d / d 2 r{H.-R-sy. 
Js 

There is also a bulk magnetic interaction 
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(48) 



The effect of stationary flow on the texture could be in- 
corporated by the dipole- velocity interaction -Fbvi but 
the flow velocities in the experiment H are so small that 
Fdv is negligible. 

The values of the many coefficients appearing above 
are discussed in Ref. [06| in some detail. However, most 
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FIG. 4. Gradient energy parameters Ag G (dashed line: 
specular wall, dash-dotted line: diffusive wall) and 2Ag2 (solid 
line) in the weak coupling approximation. The upper curves 
are for F? — and the lower ones for Ff = — 1. Fg and 
higher order parameters are assumed zero. The unit is eo£o, 
where eo = hvFN(0)kBT c . 
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FIG. 5. Surface dipole interaction parameters 62 (dashed 
lines) and 64 (solid lines) in the weak coupling approximation 
and with Ff — 0. The lower of each pair is for a diffusive wall 
and the upper for a specular wall. The unit of hi and 64 is 
5D(fc B T c ) 2 Co. 



of them have only been evaluated in the GL region so 
far. As a byproduct of our calculation of the surface 
order parameter, we can now extend the calculations of 
62, &4 and Asg to all temperatures. 

The surface gradient parameter consists of two contri- 
butions Asg = Ag G + 2AG2- The contribution Ag G is equal 
to Jy P ™/L y , the current of y-directional spin projection 
in the x direction per unit length in the y direction (L y ), 
calculated for the order parameter A^ ). In terms of spin 
current density j!P in (15) it is given by 



\ a — 
A SG — 



dzj^>(z) 



(49) 



Figure || gives the temperature dependence of both terms 
of Asg for the cases of a specular and a diffusive wall. 
Close to T c all of the parameters vanish linearly in T — T c . 
The slopes of Asg = Ag G -f-2A G 2 agree with the GL results 
of Ref. Bq] for vapor pressure; at low temperatures there 
is considerable deviation from the linear GL behavior. 
We have calculated the results for Ff = and Ff = — 1, 
the true value at vapor pressure being probably some- 
where in between. The change with Ff is rather strong 
here, since Ag G is directly related to the spin current. 

The surface-dipole coupling constants 62 and &4 are 
given by Q 

62 = t.9d / dz(A\ - 6A_lA|| + 5Af ) 
4 Jo 

bi = -9d / dz(A|, - A ± ) 2 (50) 

8 Jo 

at all temperatures. Figure @ shows these two plotted as 
a function of T/T c . Both of them go to zero proportional 



to (1 — T/Tc) 1 / 2 near T c . In this region more accurate 
values can be obtained from the GL results |3q| , with 
which these coincide. Because 62 > 264, i*sD ( M ) favors 
n perpendicular to the wall. The dependence of 62 and 
64 on F® is much weaker than that of Asg ; since the 
effect of the spin current comes to play only through the 
order parameter. For F" — —1 the values tend to be 
diminished by at most 5 percent from those shown in 
Fig. § 



B. Competing interactions and length scales 

In dealing with the pinhole model, we are mostly in- 
terested in the behavior of the texture at different mag- 
netic fields near a flat surface. In this case there are 
three relevant contributions to energy (per surface area 
L 2 ): (i) the surface-dipole energy, Fa^/L 2 ss b 2 , (ii) the 
surface-field energy, Fga/L 2 ~ dH 2 , and (iii) one related 
to the bending of the texture when there is a uniform 
perturbation at the wall (caused by i^sH, for example) 
(F G + F BK )/L 2 w ^65A G2 a J ff 2 /8 oc H. All of these have 
different field dependences, but their values turn out to 
coincide at H rs 1 mT. It is seen that at fields H < 1 
mT the constant .Fsd always dominates, and thus fi will 
be aligned perpendicular to the surface. For H > 1 mT 
the dominant surface interaction is -Fsh, and the local 
texture is determided by its minima. 

The conclusion about the relative magnitudes of -Psd 
and .Fsh is true also in general. However, to determine 
the texture at the surfaces in more complicated restricted 
geometries, the gradient and other bulk energies have to 
be minimized together with the surface energies. The 
surfaces can then be seen as perturbing the bulk texture 



10 



which would otherwise be uniform. Associated with each 
hydrostatic interaction competing with the gradient en- 
ergy, there is a characteristic length scale that describes 
the scale on which such local perturbations from unifor- 
mity will decay, or heal p7p7| ] . The stronger the compet- 
ing interaction, the shorter is the corresponding healing 
length. Comparing these healing lengths with the spatial 
scale I 3> Co of the container gives qualitative information 
on the form of the texture. 

Most importantly, in the presence of a magnetic field 
Fdh and Fq define a length £h = \/65Ag2/ (8aH 2 ) ex 
H^ 1 J36|. Comparison of Fd and Fq gives another one, 
the dipole length £ D = V^WffD^ w 10 /mi > £ , 
related to possible perturbations of the rotation angle 
9 from its equilibrium value 9q. As stated above, we 
assume there to be no interactions present which could 
force such perturbations, and therefore £d plays no role. 
In any case, for most practical purposes £h, I 3> £d, and 
such a perturbation would decay fast on the scale of /. 

For very small magnetic fields £h ^S> I and £h is thus 
also not a relevant length scale. The only important hy- 
drostatic interactions in this case are the surface-dipole 
energy Fsd and the gradient energies Fq and Fgc- Now 
the interesting question is simply whether the walls of the 
container are separated by long enough distances for Fsd 
to be essentially minimized, or small enough distances 
for no textural variation to occur at all — the minimum 
configuration of -Fbtot = Fq + -Fsg- An elementary esti- 
mation of the length scale I at which there is a transition 
from one behavior to another gives I ex AG2A&2 — ^4), 
where the constant of proportionality is of order unity. 
Since 64 < &2, we can drop it and define a surface-dipole 
length £sd = Ag2/^2- Using the numerical values calcu- 
lated above we get at zero pressure 
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VIII. ISOTEXTURAL JOSEPHSON EFFECT 



The Josephson energy ra) of a pinhole depends nontriv- 
ially on three parameters: the phase <fi and two parame- 
ters describing ipy = R^R^. In addition there is depen- 
dence on the surface scattering (diffusive vs. specular), 
on the temperature T/T c and on W/D. Below we can 
present only some representative parts of this parameter 
space. In this section we plot isotextural current-phase 
relationships, where ipij is assumed constant in J(<j>, ipg). 
The possible <f> dependence oiipij is considered in Sec. [LX| . 

The pinhole coupling has the maximal symmetry (|6J) 
independently of the shape of the hole as long as W/D = 
0. In terms of n this means, for example, Fj{(j>, n L , n^) = 



i.R *l\ __ f 3 (<I), -n L ,-n R ) 
mass current. The mass current is given in units of 



Jo = 2rri3V-FN(0)kBT c A o , where A is the total open 
area of one or more pinholes. In most cases only the 
phase differences in the range [0, n] are shown due to the 
symmetry J(2ir — <j>) = —J((j>). All plots are made for 
Ff = 0. Tests with Ff = —1 show no qualitative differ- 
ences and at most a few percent quantitative difference 
in J{4>) even at the lowest temperatures. 



A. Spin-rotation axes perpendicular to wall 

We first study the case where the spin-orbit rotation 
axes h L:R are perpendicular to the intervening thin wall. 
This situation is realistic if the external magnetic field is 
small enough (H < 1 mT) and if there are no other walls 
with different orientations nearby. Four different fi con- 
figurations are then possible, namely the combinations 



n 



L,R _ 



±z, where z is normal to the wall. These give rise 
to two different current-phase relations corresponding to 



parallel (n = fi ) or antiparallel (n 



n 



= ±z) 



Fj((/>, n K , n- L ) = Fj(<j), — n , — n ), and similarly for the 



situations. The parallel case is actually more general, be- 
cause the fi vectors need not be perpendicular to the wall 
to still give the same J(4>). The current-phase relations 
are shown in Figs, mq corresponding to three different 
surface models. 

First we consider the case of constant order parameter. 
Here A(k, r) is not calculated self-consistently using any 
boundary condition, but instead is assumed to have its 
constant bulk form Ak all the way to the wall (no pair 
breaking). This is the case discussed by Yip pi], and 
the current-phase relations shown in Fig. rc| are exactly 
the same as obtained by him. The parallel case is well 
known: the same result was first obtained by Kulik and 
Omel'yanchuk for microbridges in s-wave superconduc- 
tors, and it was subsequently generalized to the case of 
3 He by Kurkijarvi |16J . The new feature found by Yip 
is seen in the antiparallel case: very close to T c the J((f>) 
is sinusoidal, but at temperatures below about 0.5T C the 
7r state appears: a strong kink and a new zero of J((/>) 
appear around cj) = tt. 

The self-consistent surface models lead to considerable 
suppression of the order parameter at the wall. As a 
consequence, the J(^>)'s are different for a specular sur- 
face and for a diffusive surface, shown in Figs. M and 
H. Both surface models result in qualitatively similar 
J(0)'s. The parallel current-phase relations look simi- 
lar as in Yip's case, although their critical currents are 
slightly reduced. A clear difference is seen in the antipar- 
allel cases: Firstly, the whole </(</>) appears to be shifted 
by 7r so that on the phase interval [0, it] the current is 
mostly negative. Secondly, the J(<f>) remains sinusoidal 
down to very low temperatures. An additional kink be- 
gins to form only at around 0.2T C . Now the kink is also 
shifted from (j> — tt to <p — 0. We continue to call it a 7r 
state because it represents a local minimum of Fj (</>) that 
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FIG. 6. Isotextural current-phase relations for a pinhole in 
a wall with a constant order parameter on both sides. The top 



panel corresponds to n L 

W/D = 0. 



n and lower to n 



\R 



±z. 



is shifted from the from the global minimum of Fj (0) by 
the phase difference 71". 

Figure shows the critical currents J c and the possible 
additional extrema of </(</>) as a function of temperature. 
For parallel fi vectors such a plot has been published in 
Rcf. p4[ . but there the pinhole results were erroneous. 
Close to T c , J C (T) oc 1 — T/T c for a constant order pa- 
rameter and a specular surface, and for a diffusive surface 
J C (T) oc (1 — T/T c ) 2 , as expected from previous calcu- 
lations |23|. The critical current for a constant order 
parameter is always the highest and for a diffusive wall 
the lowest. For antiparallel n vectors the roles change: 
the constant order parameter case has the lowest J c , due 
to the strong cancellation between different quasiparti- 
cle directions, but the negative extremum around <p = 7r 
is nearly as pronounced as the positive one. It is clearly 
visible that the other extrema appear only at much lower 
temperatures for diffusive and specular surfaces. 

The dotted lines correspond to the high-temperature 
approximations obtained from Eqs. (|4l|) for a diffusive 
wall, i.e., (2m3 / 'H)(a + 2/3) for parallel and {2m^/K)\a — 
(7/4)/3] for antiparallel n's. These lines follow the correct 
critical currents very well down to temperatures around 
T = 0.4T C . The current-phase relations show some de- 
viation from the sinusoidal form at temperatures above 
0.4T C in specular and diffusive cases, but the deviation is 
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FIG. 7. Isotextural current-phase relations for a pinhole in 
a specularly scattering wall. The top panel corresponds to 

-n fl = ±z. Ff = W/D = 0. 



much smaller than for a constant order parameter. 



B. Other orientations of spin rotation axes 

In order to study other orientations of n L < R we con- 
sider first the case of relatively large magnetic fields, 
H > 1 mT. The configurations, which minimize the sur- 
face magnetic interaction Fsu (J47]), depend on the an- 
gle 6h between H and the wall normal s = z. Our 
results, obtained using the self-consistent order param- 
eter for a specular and a diffusive wall, are qualita- 
tively similar to Ref. p0|| , but differ in details. Figure 
[l0| shows the four different </(</>) ' s which are possible at 
Oh = 0.45-7T, corresponding to four different n L ' R config- 
urations A A, AB,AC and CD, as defined in Ref. JL0|. 
Also shown is the dependence of J {(f) on 8h, when the 
n L ' R are in the AB configuration. For Oh = the AB 
configuration gives the antiparallel case studied above. 
These should be compared with Figs. 4 and 5 of Ref. 
|lQ|, where the same plots were given for the constant 
order parameter. It can be seen that it states are not 
uncommon at low temperature T ~ 0.1T C . 

For a systematic study it seems to be more economic 
to specify n L ' R directly. Let us consider the case where 
the polar angle of n L is changed but n^ is kept con- 
stant: h L ■ z = cos7/ L and n^ = z. It can be seen in 
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FIG. 8. Isotextural current-phase relations for a pinhole 
in a diffusely scattering wall. The top panel corresponds to 
n L = n R and lower to n L = ~n R = ±z. Ff = W/D = 0. 



Fig. hip that J{4>) differs essentially from the sin (j> shape 
around angles i] L sw 0.467T. The current is nearly pro- 
portional to sin(20) at ?y L = 0.467T, and it states are 
present in the range r\ L « 0A2tt . . .0.50ir at tempera- 
ture 0.4T C . The top and bottom solid lines in Fig. Ob 
show the extrema of J{4>) as a function of r\ L . The 
middle solid line shows the prediction of the tunneling 
model. The tunneling model (f7|) has always sinusoidal 
isotextural J(<f>) = J c sin</>, where J c = (2m3/h)E c and 
E c = aR^R^+PiR^R^+R^R^). Therefore, it has 
only a single extremum value J c (in the range < <f> < 7r). 
Fig. [ll| b shows that the n states (where two extrema 
appear) take place around the configuration where J c 
changes sign by going through zero. 

Fig. [ll]b shows the extrema at three different temper- 
atures. At high temperatures the range where n states 
occur is very narrow. Also, in configurations showing a 7r 
state, both the negative and positive extrema of J(4>) are 
much reduced relative to the maximal critical currents 
shown in Fig. ||a. Outside of the range where E c rs 0, 
</(</>) is nearly sinusoidal and the critical current is well 
predicted by the tunneling model. With decreasing tem- 
perature the range of the configurations showing a it state 
widens and the extrema get larger. These results are not 
restricted to the case where n^ = z but are valid for 
general configurations of h L,R . This can be seen as a 
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FIG. 9. Critical currents as functions of the temperature. 
The lines (with and without markers) are theoretical pinhole 
results and separate points (x, *) are experimental results M 
with Jo = 213 ng/s (Sec. tX[). For theoretical results the upper 
and lower panels correspond to parallel and antiparallel n vec- 
tors, respectively. For experimental results they correspond 
to the H and L states. The signs (+, — ) denote the negative 
and positive extrema of J{4>) that appear in the case of an- 
tiparallel n's. The experiment shows two extrema in both H 
and L cases, whose signs are unknown. The calculations are 
for Ff = W/D = 0. 



manifestation of a general difference between tunneling 
junctions and weak links, as discussed in Ref. |}8|. In sit- 
uations where a tunneling supercurrent is prohibited by 
symmetry, there can still be a small supercurrent flow- 
ing in a corresponding weak link, although with some 
restrictions on the form of </(</>). 



C. Nonzero W/D 

Let us consider quasiparticle scattering inside a pin- 
hole caused by a finite aspect ratio W/D. In a p-wave 
superfluid this can lead to a current-phase relationship 
which is more complicated than some average of the ones 
obtained using tunneling model and direct transmission. 
Consider a deflected trajectory (Fig. 0) which is trans- 
mitted but the momentum direction is nearly reversed. 
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FIG. 10. Examples of current-phase relations in high mag- 
netic field, (a) The four possible n configurations which result 
in different current-phase relations for magnetic field angle 
Or = 0.457T. The nomenclature follows the definitions of Ref. 
[Il0|. (b) The effect of varying 0h in the AB configuration. 
The values of Oh are given in the legend. The curves are for a 
pinhole in a diffusive wall at T = 0.1T C with Ff = W/D = 0. 



Assume n L = ri R . If the phase difference 4> is zero, the 
quasiparticles of 3 He see effectively a change in the sign of 
the order parameter A(k) along such a trajectory. How- 
ever, if (j> is near tv, such quasiparticles see effectively a 
constant order parameter. This means that the scatter- 
ing reduces the energy at <fi — n relative to the energy at 
<j) = 0. This could work as a mechanism for the formation 
of 7r states even when the spin-rotation axes on the two 
sides are equal. This mechanism is closely related to the 
7r state mechanism of Ref. llG] and the effects discussed 
in Ref. |39] for sd contacts. 

The effect of the scattering on the current-phase rela- 
tion is shown in Fig. [12| The main change is the decrease 
of the critical current. Contrary to our expectation of 
forming a 7r state, a dip develops at <f> w 7r/4 with in- 
creasing W/D. However, this happens only in the region 
of W/D where the distribution ( gg ) has probably ceased 
to be valid. For small W/D the only effect is to reduce 
the critical current, and the relative initial decrease does 
not seem to depend essentially on the temperature. The 
results for specular scattering are very similar. 
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FIG. 11. Change of isotextural current-phase relation 
and extremal currents (b) as a function of r\ L 
for £i R = z. In (a) the current-phase relations are shown 
at intervals Ar/ = 0.027T at temperature 0.4T C . In (b) the 
extrema of J{4>) (top and bottom curves of each type) are 
compared with the tunneling model results (middle curves) 
at temperatures T/T c = 0.1, 0.4, and 0.6. The angles r/ 
and ty — tjq corres pond to possible bistable states, as will be 
discussed in Sec. XA. The figure is calculated for diffuse 
scattering and F" = W/D — 0. In the specular case the 
relative deviation from the tunneling model is slightly larger 
(-20%). 
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FIG. 12. Current-phase relations for diffusive scattering in 
the aperture for T/T c = 0.1, Ff = and W/D = 0.0 . . . 0.7 
in order of decreasing critical current. The n vectors on are 
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IX. ANISOTEXTURAL JOSEPHSON EFFECT 

In the previous section we assumed that the n texture 
remains constant when the phase difference </> is changed. 
In this section we study the anisotextural Josephson ef- 
fect where the texture is allowed to change as a function 
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of (j) [^lj . We demonstrate the anisotextural effect using 
a simple model. The possible existence of either isotex- 
tural or anisotextural Josephson effects in the Berkeley 
experiment J5| is discussed in Sec. [XJ 



A. Phenomenological model 

We consider a planar wall separating two half-spaces 
of 3 He-B. In the absence of perturbations, the surface in- 
teraction (El) fixes a constant texture on both sides. The 
orientations on left and right hand sides are denoted by 
n^ and n^,, respectively. They both are either parallel 
or antiparallel to the normal z of the wall: n^ = ±z 
and n^ = ±z. Let us now place a weak link in the wall. 
The Josephson coupling energy Fj (|3J) may now favor a 
different orientation h L ' R ^ ±z at the junction. Such a 
change is opposed by a "rigidity energy" , which consists 
of gradient energy and possibly other textural energies 
([43|)-(ff8|). We model the rigidity energy by a quadratic 
form 



F Tis = J L (v L -V^) 2 +l R (v R -V^ 



(52) 



Here r) L ' R and rj^ R are the polar angles of h L ' R and 
n^ R , respectively. In the case of symmetric left and 
right sides, the stiffness (or rigidity) parameters j L > R are 



7), but the more general form will be 



equal (j L = 7 J 
useful later. 

The texture is now found by minimizing the free energy 



F = F, 



,n 



,n«) + F rig (n-,n«,n^,n£) 



(53) 



with respect to n L and n^. If the stiffness parameters 
jL,r are su fg c i en tly small, this will lead to a dependent 
texture. This is easiest to see using the tunneling model 
(0) for Fj. If cos0 > 0, the minimum Fj = — (a + 
2/3) cos<p is achieved by n L = h R . In the opposite case 
cos <fi < the minimum Fj = — (a — 2/3) cos <fi (assuming 
a < /?) is achieved by h L > R = (=px + y =p y/3z)/y/5, for 
example. Neglecting the rigidity (|5^), this leads to a 7r 
state with a piecewise sinusoidal current-phase relation 



J(4>) = 



(21713/ h) [oi + 2/3) sin 4* for cos > 
(21713/ h) (a — 2/3) sin for cos < 



(54) 



This ideal current-phase relation is smoothed by finite 
ryL,R^ wif n increasing stiffness, the texture changes less 
as a function of cf>, and for j L > R ^> a, (3 the current-phase 
relation reduces to the isotcxtural one. 



B. Model parameters 

Instead of using the tunneling model, we now assume 
the weak link to consist of an array of pinholes. The 
phase c/> an d n L ' R are assumed to be constants over the 



array. Consequently, we can use the single-pinhole re- 
sults for current ( p8| ) and coupling energy (Ejj) simply by 
replacing the open area A by the total area of the array 
A times k, the fraction of area occupied by the holes. 

We start to estimate 7 including first only the gradient 
energy: 



F ■ 

-frig 



o A G2 



d 3 r I6(a i n J )(a i n i ) 



-(\/3V-n + VZh- V x hf 



(55) 



This can be obtained from the sum of gradient energies 
(|46|) and (EJ) assuming A S g = 4A G 2 and A G i = 2A G2 ; 
see also Ref . [J37J . We shall assume that n varies only in 
one plane, the xz plane, for example. This allows one to 
describe n by its polar angle r\ alone: n = cos rji + sin r/Sc. 
In addition, we can assume r) to depend only on the radial 
direction r from the center of the aperture array. With 
these assumptions the energy (|55|) on one side simplifies 
to 



^ S = ^A G2 ldrr\d r r,f 



(56) 



This is minimized by a function of the form r\ (r) = A/r + 
C, but to avoid a divergence at r — 0, we have to cut off 
the integration at some r, for example at the radius of 
the array R = \J 'A/n. As a result one finds the form (|52 
with the stiffness parameter 



7 



50tt 



A G 2-R- 



(57) 



The rigidity has also a contribution form the surface- 
dipole interaction ( [44] ) . It can be shown that this correc- 
tion to 7 is proportional to R/S,sb when the radius of the 
array R is small compared to £sd (pl|)- We assume this 
to be the case, and therefore neglect the small correction. 

In order to get numerical values for 7 ( p7| ) we use the 
Berkeley array, where R ~ 0.11 mm pi. The tempera- 
ture dependence is given by A G 2, shown in Fig. 0. The 
resulting j(T) is plotted in Fig. |[ We use this value as a 
standard, to which the parameters j L - R used in our cal- 
culations are referred to. Note that the reference value 7 
of the rigidity energy is larger by one order of magnitude 
than a and (3, which give the magnitude of the Josephson 
coupling energy (Fig. 0). 

In order to calculate the current, the free energy 
( p3| ) should be minimized with respect to three angles 
parametrizing r\. L,R : rj L ,ii R and the relative azimuthal 
angle \- The results of such a minimization for j L — 
~/ R = O.I7 are shown in Figure 111. In creating these 
curves, we proceeded from left to right, using the mini- 
mum angles (r] L , rj R , x) of each cf> step as the initial guess 
for the next step. The panels on the left are for the par- 
allel case (77^ = "Hx — 0)- For (j) from to approximately 
tt/2, the vectors h LR remain exactly perpendicular to 
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FIG. 13. Anisotextural current-phase J(0) and 



en- 
ergy-phase ^(0) relations for bistable states 77^ = (a,c) and 
r]^ — tv (b,d). The curves correspond to different tempera- 
tures at intervals of 0.1T C . The arrows indicate the fi orien- 
tations of in the zero branch, where n L — n^ and & R = n^. 
The rigidities on both sides are equal, 7 L = 7^ = O.I7 and 
?7oo — 0. The weak link parameters and 7 (Fig. H) are eval- 
uated for the Berkeley array (Sec. Kf). The wall is assumed 
diffusive, Fi° = W/D = 0. 



the wall. At low temperatures, a discontinuous jump to 
another branch of J((f>) occurs at around <f> = tt/2, where 
n LR are tilted form their original perpendicular posi- 
tions. With increasing (j) there is a jump back from this 
7r branch so that for <f) > 3tt/2 the minimum solution 
again corresponds to perpendicular ii's. The panels on 
the right are for the antiparallel case (77^ = tt, 77^ = 0). 
Here the tt branch, where n's are tilted, occurs not at 
<f> rs tt but at <j) w 0. The fact that the curves are not 
(anti)symmetric with respect to <j> — tt indicates that the 
jumps between different branches are hysteretic. 



We have studied the effect of increasing 7 1 



7 



from 



the value O.I7. We find that the tt state first disap- 
pears in the antiparallel case and then also in the parallel 
case. For example, the tt state in the parallel case [de- 
fined as positive J'(tt)] disappears when j L < R rs O.27 at 
T = 0.4T C . The anisotextural effect on the current-phase 
relation still continues up to j L ' R sa O.77. Since the cou- 
pling ( (40| ) scales with A = An = ttR 2 k and stiffness ( |57| ) 
with R, we get the following necessary condition for the 
appearance of tt states in the parallel case at tempera- 
tures T > 0.4T r 



Rk > 0.5 /im. 



(58) 



Thus, the larger the radius R of the array and the higher 
the ratio of the open area k, the better are the chances of 
realizing the anisotextural tt state. At higher tempera- 
tures it will be more difficult, due to the relative strength 
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FIG. 14. Anisotextural J(4>) for different asymptotic an- 
gles r/^ at temperatures 0.2T C (a) and 0.4T C (b). The other 
parameters are 77^ = 0, 7 L = O.37, and 7^ = 7. The dif- 
ferent curves correspond to r?^, at intervals of 0.1-7T. Other 
parameters are the same as in Fig. [L3J. 



It is interesting to compare the condition (pq) with the 
Ginzburg-Landau calculation in a single large hole [pTj. 
There tt states could be seen for hole radii R > 5.5£gl- 
Setting k = 1 and extrapolating £gl(0.4T c ) = 8.8 nm 
fl36| , we have the condition Rk > 0.48 /jm in surprising 
agreement with the pinhole result above. 



C. Asymmetric case 

We go slightly beyond the simple model of an infinite 
planar wall introduced in Sec. IXA. Firstly, we allow 
the asymptotic directions n^^ to be arbitrary. Secondly, 
the stiffness coefficients j L > R can be different. In Fig. [IJ 
we have studied different values of 77^ while 77^ = 0, 



7 = O.37, and 7 = 7. Here the extreme curves repre- 
sent parallel and antiparallel n^ R . It can be seen that 
neither of them have a tt state at 0.4T C and only the 
parallel case has one at 0.2T C . However, at intermedi- 
ate 77^ the tt states still persist in a wide range of 77^. 
The current-phase relationships are not hysteretic at high 
temperature, but hysteresis develops at lower tempera- 
tures. Figure O should be compared with the corre- 
sponding isotextural Fig. O. In the isotextural case the 
tt states occur only in limited range, and there is no hys- 
teresis. 

The states with 77^. = tt/2 ± constant are expected to 
be degenerate in the absence of the Josephson coupling 
(t7oo = 0). The currents and energies for one pair of 
such bistable states are shown in Fig. [Tjj The impor- 
tant difference to Fig. [13] is that the curves are smooth 
and there is no hysteresis. In the case of parallel and an- 
tiparallel fine's the symmetry is spontaneously broken in 
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FIG. 15. Anisotextural J((j>) and F(cj>) for bistable states 
ry^ = 0.3-7T (a,c) and r^ = 0.7ir (b,d). The textures have 
different rigidities r y L = 0.37, J R = 7 an d ?7oc = 0. Other 
parameters are the same as in Fig. O. 



the 7r branch, whereas the tilted n^ already breaks the 
symmetry, and thus the it state can develop continuously. 



ference between the infinities can be obtained by adding 
(2m,3 / h) J / (ir Rp s ) to 4>, where the correction assumes a 
radial flow outside of the radius R. Correspondingly the 
total energy gets an additional contribution J 2 /(2irRp s ). 
These rescalings do not affect the results qualitatively. 

The anisotextural model, as described above, assumes 
a vanishingly small external magnetic field. Qualita- 
tively, it is easy to see what the effect of a strong magnetic 
field would be. In a field H ^ 1 mT, the strongest inter- 
action affecting the texture is the surface-magnetic term 
Fsh (f§7])- If the coupling energy scale is much smaller 
than this, then the texture will be fixed to some mini- 
mum of -Fsh and will not depend on the phase difference. 
A strong magnetic field therefore suppresses any aniso- 
textural 7T state and only the isotextural mechanism re- 
mains. 

In order to estimate the critical field, we equate the 
change in magnetic surface energy AF§h ~ dAH 2 at 
the junction with the gain of energy AE between the tt 
state and "0 branch" at = n. Thus H c « yjAE/(dA). 
For the Berkeley array H AE can be estimated from 
the energy-phase graphs of Figs. [l3| and |l5] or from the 
experiments. At T rs 0.45T c the values are on the order of 
AE w 0.1 . . . 0.5 aJ, which yields the order of magnitude 
H c rs 10 . . . 50 mT. 



D. Discussion 



The results of Figs. [13-h5l contain both the isotextu- 
ral and anisotextural mechanisms of n states. However, 
practically the same results can be obtained all the way 
down to T = 0.4T C by using the tunneling model Fj 
(R), which excludes the isotextural n state. As discussed 
above, the tunneling model fails at high temperatures 
only if the Josephson energy is close to zero. The mini- 
mization procedure seems to avoid such a situation, and 
thus the tunneling model gives a good description of 
the anisotextural pinhole array at temperatures above 
T w 0.4T C . 

Here we comment our earlier calculations of the aniso- 
textural Josephson effect in Ref. fllj| . Firstly, there the 
tunneling model was used instead of the general pinhole 
result. Secondly, the scattering within the hole was par- 
tially taken into account by reducing the transmission 
by factor p(i?) (Pjj), but the deflected trajectories were 
completely neglected. Thirdly, the coupling parameters 
a and were scaled by different factors, which is not 
possible when using the general pinhole result. Taken 
together, these explain why Fig. O is different from Fig. 
2 in Ref. @. 

To be accurate, we should add to the preceding cal- 
culation the effect of the flow in the macroscopic region 
far away from the weak link. Above the phase differ- 
ence 4> is defined between the macroscopic-mesoscopic 
borders of the two sides (Fig. [j]). The "true" phase dif- 



X. ANALYSIS OF THE BERKELEY 
EXPERIMENT 

We now turn to an analysis of the Berkeley experi- 
ment H . There the weak link consists of a square array 
of 65 x 65 holes. They were etched in a 50 nm thick sil- 
icon chip with a hole spacing of 3 /iin, making the area 
of the array 195 /im x 195 /im. The holes were nom- 
inally squares 100 nm x 100 nm. However, flow resis- 
tance measurements in the normal state seem to indicate 
a somewhat larger apertures 115 nm x 115 nm |4fJ], and 
these larger values are used in all numerical estimates in 
this paper. A sketch of the experimental cell is shown in 
Fig. [16. The scattering properties of 3 He quasiparticles 
from the silicon chip are not known, but most surfaces 
are generally believed to be diffusively scattering. The 
magnetic field is believed to be small, H <C 1 mT, and 
the pressure is bar. We further assume the system to 
be in thermal equilibrium. 

The central experimental findings are the bistability 
and the existence of tt states || . Bistability means that 
the system can randomly choose between two alternative 
states, characterized by high (H) and low (L) critical cur- 
rents. Both of these states show tt states. The measured 
extremal currents are plotted in Fig. 0. 

There are two major difficulties in applying the the- 
ory presented in this paper to the Berkeley experiment. 
Firstly, since £q — 77 nm, the holes are too large to be 
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weak link 




FIG. 16. Part of the experimental cell in Ref. ||. Most of 
the structure has rotation symmetry around the vertical axis. 
The "pillbox" (a) between two flexible membranes (diameter 
1.27 cm, distance h — 0.14 mm) forms a volume that is con- 
nected to the rest of 3 He volume only through the weak link. 
The weak link is made in a silicon chip attached to the lower 
membrane. The chip is /12 = 0.5 mm thick and has a square 
window opening from 250 fim x 250 /im to some 0.95 mm x 
0.95 mm at the lower chip surface. The figure is based on a 
drawing supplied by S. Pereverzev. 



pinholes. Also, the holes are too small for the Ginzburg- 
Landau calculations to be reliable pLO]. We use the 
pinhole model because more accurate calculations would 
be much more demanding. Moreover, due to the approx- 
imate nature of our pinhole calculations for finite aspect 
ratio W/D, we will only use the pinhole theory in the 
limit W/D = 0, although experimentally W/D fa 0.4. 
Another reason for using W/D = is that the measured 
critical currents are clearly larger than calculated for pin- 
holes with W/D = 0.4 (Sec. |VIIIC| ). The second major 
difficulty is that the cell is complicated, and its dimen- 
sions are on the same order of magnitude as £sd ( pl|) ■ In- 
stead of a proper calculation of the texture, we will make 
some simple estimates and introduce one adjustable pa- 
rameter. 



A. Isotextural Josephson effect 

Here we consider the case where the Josephson cou- 
pling can be considered as a weak perturbation, which 
does not affect the texture on either side of the weak 
link. For small magnetic field, the dominant orienting 
effect on n comes from the surface-dipole energy (H). In 
region (a) of Fig. Wty this clearly favors the uniform tex- 
ture n = ±z, where z is along the axis of the cell. The 
situation is more complicated on the other side of the 
junction. The axial orientation is preferred in the wide 
cylindrical region (d). The narrow cylindrical region (e) 
favors n _L z. The tendencies from regions (d) and (e) 
compete in region (c) below the chip, and affect the tex- 
ture in the window region (b). Let us assume that the 



minimization of textural energies (J43|)-([48]) favors at the 
junction a particular orientation of n = rip with po- 
lar angle rtf. Then there must be a degenerate textural 
state with h L — ~Uq, which corresponds to the polar 
angle r\ L = tt — t]q. There can be additional degeneracy 
with respect to the azimuthal angle. 

When we add the pinhole Josephson coupling, the de- 
generate states above give rise to precisely two different 
current-phase relationships if t]q ^ tt/2. These two are 
obtained by studying configurations with r\ = t]q and 
i] L = tt — t]q with a fixed n^ = z. Configurations of this 
kind were studied above in Fig. Flh. They can in prin- 
ciple explain both the bistability and the existence of tt 
states. Quantitative comparison with experiment gives, 
however, very poor agreement. The problem is to fit the 
four experimental points in Fig. |E)| at any temperature 
with the four extrema in Fig. |ll| b at angles rj^ and tt — t]q 
using ?7q as a fitting parameter. One such construction 
at T — 0.4T C is shown by vertical lines in Fig. [ll|b. In 
that case, one finds the L state with a tt state, but then 
no 7r state appears in the H state, and all currents are 
more than by a factor of 2 too small. Alternatively, if one 
tries to fit the critical current in the H state, one has to 
approach the parallel and antiparallel states, where no 7r 
states appear at this temperature (Figs. Rfl). 



B. Anisotextural Josephson effect 

For the anisotextural Josephson effect we have to cal- 
culate the rigidity energy (p2|). The region (a) is much 
thinner than £sd and therefore we have to consider a 
two-dimensional texture instead of the 3D in Sec. LXB| . 
In the present case the surface interaction is important. 
We find that -f R is proportional to Ag2^/ ln(V 7i£sd / R) 
instead of Ag2-R (pd), but these two happen to be of the 
same order of magnitude. Thus the texture in region (a) 
is rather stiff, and we assume 7^ « 7. This makes the 
texture essentially fixed. 

In order to analyze the other side (L), we first consider 
a conical region between two radii, R\ and R 2 . Oth- 
erwise, we use the same approximations as for the half 
space in Sec. KB We get the rigidity energy 
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R1R2 
Ri — R\ 



[r,(R 1 )-r,(R 2 )} 2 , 
(59) 



where 9 is the opening angle of the cone. This reduces to 
the previous result ( p7| ) in the limit 9 = tt/2 and R2 — 
00. We now consider the regions (b) and (c) as three 
consecutive conical regions. The middle one corresponds 
to the conical region of the chip with 9 — 9 W fa 38°, 
the two others have 9 — tt/2. We neglect the effect of 
regions (d) and (e) and assume that region (c) extends 
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to infinity. For the combined stiffness of these regions we 
get the estimate 
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(60) 



where R\ and i?2 are the radii of the inner and outer 
edges of the conical region measured from the weak link. 
Substituting the numerical values we find j L m O.3I7, 
where 7 is the half-space value (pTJ). Thus this side is 
considerably softer than the other. 

The calculations for t he p resent parameters have al- 
ready been done in Sec. IXC . The results are presented 
in Fig. [lj for different angles rj^. A representative pair 
of bistable states is presented in Fig. [l5|. The calcu- 
lated current-phase relations are very similar to those 
found experimentally g. The 7r states are present in 
both bistable states. The critical current in the H state 
(identified with r]^ w 0.37r) is very close. Some differ- 
ences can be resolved. For example, the 7r state is too 
strong in the theoretical H state, and the critical current 
in the theoretical L state is slightly too large. The only 
fitting parameter here is the textural angle ij^, and it 
can be seen from Fig. |l4] that rj^ = (0.5 ± 0.2)7r roughly 
represents the best overall fit to the experiments J5| . Tak- 
ing into account that the experimental apertures are not 
pinholes and that only one fitting parameter is involved, 
the agreement between the anisotextural theory and the 
experiment is amazingly good. 

The anisotextural Josephson model gives several pre- 
dictions that can be tested experimentally. The shape of 
the current-phase relationship crucially depends on the 
number of parallel apertures and also on their spacing. It 
depends on the geometry of the cell surrounding the weak 
link. It depends on the magnetic field. The current-phase 
relationships become hysteretic at low temperatures. All 
these dependences are either absent or very different in 
the isotextural model. The dependences can be quanti- 
tatively extracted from the theory presented above. 

None of the predictions have been studied experimen- 
tally yet, excluding possibly the hysteresis. The first 
paper j| reported discontinuous jump to the it branch 
at low temperature 0.28T C . The later paper || reported 
continuous current-phase relationships but only at higher 
temperatures > 0.45T C . However, there was a change 
in the experimental set up between these observations, 
which may have affected the results. The hysteresis 
should also show up as additional dissipation, but no 
detailed theory yet exists Elf . 



XI. CONCLUSIONS 

We have presented above a fairly complete study of the 
dc Josephson effect in 3 He-B using the pinhole model. 



We have derived a general energy functional for the pin- 
hole coupling energy. A computer program has been con- 
structed to calculate the coupling energy and currents in 
a general case. We have plotted the isotextural current- 
phase relationships for various cases. Besides the mass 
current there is also a spin-current, but that has not been 
examined in this work. In contrast to most pinhole calcu- 
lations, we consider also a finite aspect ratio of the hole, 
but no extensive studies have been made because of the 
approximate nature of the model. We have calculated 
some surface parameters of 3 He-B. 

We have studied the anisotextural Josephson effect. 
The previous tunneling model calculations are general- 
ized to arrays of pinholes. The general conditions for 
the anisotextural effect are studied. It is found that the 
anisotextural Josephson effect depends sensitively on pa- 
rameters like the dimensions and the number of holes, 
the surrounding geometry, and the magnetic field. 

The theory is applied to explain the experimental ob- 
servations made at Berkeley. We compare the exper- 
iment with both iso- and anisotextural models. Both 
mechanisms can in principle explain the bistability and 
the 7r states. In quantitative comparisons there is one 
adjustable parameter describing the texture. Quantita- 
tive comparison with the isotextural model gives poor 
agreement, but a very good agreement is found with the 
anisotextural model. New experiments should be made 
to confirm the identification of the anisotextural Joseph- 
son effect. 

Experiments on the Josephson effect in 3 He-B have 
also been done using a single aperture H. Both n states 
and multistability are observed. The pinhole theory can 
hardly be applied to this case because the aperture is 
much larger than the coherence length £o- Also, the 
distinction between iso- and anisotextural effects is not 
well defined for a single aperture. The Ginzburg-Landau 
calculations should be more accurate here, but unfortu- 
nately they have been done systematically only for par- 
allel n's |j. 
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